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Abstract 

An identity that relates multipolar solutions of the Einstein equa- 
tions to Newtonian potentials of bars with linear densities proportional 
to Legendre polynomials is used to construct analytical potential-density 
pairs of infinitesimally thin bars with a given linear density profile. By 
means of a suitable transformation, softened bars that are free of singu- 
larities are also obtained. As an application we study the equilibrium 
points and stability for the motion of test particles in the gravitational 
field for three models of rotating bars. 
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1 Introduction 

Bars are a common self-gravitating structure present in disc galaxies. About 
50% of such galaxies are stongly or weakly barred, including our Milky Way 
[U |2]. The effect of a weak bar is usually represented as a potential in 
cylindrical coordinates in the form (p) = cos(2</?) [3]. In the case 
of strong bars, the only exact, self-consistent models known are those of 
Freeman [JJ, although they present some unrealistic features for bars. In 
studies of orbits involving strong bars, they are often modelled as homo- 
geneous ellipsoids [SI [6] or inhomogeneous prolate spheroids [TJ EJ [U HO) . 
These mass distributions have a finite extent. Long & Murali [11] discuss 
a simple method to generate analytical potential-density pairs for barred 
systems that extend to all space. 
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In this paper we construct analytical potential-density pairs for infinites- 
imally thin and 'softened' bars that can be expressed solely in terms of 
elementary functions. The starting point is an identity that relates multi- 
polar solutions of the Einstein equations to Newtonian potentials of bars 
with densities proportional to Legendre polynomials. These bars can then 
be superposed to generate other bars with a desired density profile. We 
also use the method of [TT] to soften the infinitesimally thin bars. This is 
presented in Section [2j In Section [3] the potentials for barred systems are 
used to study an aspect of the motion of test particles in uniform rotating 
bars, namely the equilibrium points (Lagrange points) and their stability. 
We relate the properties of the equilibrium points to the mass distribution 
of the bar models. Section U] is devoted to the discussion of the results. 



2 Bars with variable densities 

In this section an identity derived by Letelier [12J will be used as a starting 
point to construct potential-density pairs of bars with various linear density 
profiles. The Newtonian potential of a bar of length 2L with linear density 
X(z) located symmetrically along the z-axis is 

* = -c[ L , Mz ' )dz ' . (1) 

J ~ L ^R 2 + (z- z'f 

where G is the gravitational constant. Letelier [12] found the following 
identity: 

n (t\P ( \ 1 f L p n( z '/ L ) dz ' M 

Qn{0 P n\V) = 7J / i =, I 2 ) 

1 J - L ^R 2 + (z- z'f 

where P n and Q n are, respectively, the Legendre polynomials and the Leg- 
endre functions of the second kind and (£, rj) are the spheroidal coordinates 
related to the cylindrical coordinates (R, z) through 

£ = (R 1 + R 2 )/(2L), V = (R 1 -R 2 )/(2L), (3) 



R 1 = ^/R 2 + (z + Ly, R 2 = y R 2 + (z — L) , (4) 

with £ > 1 and — 1 < rj < 1. Comparing equations ([I]) and ([2]), and intro- 
ducing the mass M, we see that relation ([2]) represents a family of bars with 
linear density 

M 

K(z) = ^Pn(z/L), (5) 
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associated with a potential $ n = —GMQ n (£)P n (rj) / L. 

Since the Legendre polynomials form a complete set of functions, the 
members of the family ([5]) can be superposed to generate potential-density 
pairs for bars with a prescribed density distribution. The simplest case is 
the bar with constant density 



A - M 
A °-2L' 

whose potential can be expressed in cylindrical coordinates as 
GM 



(6) 



2L 



In 



z-L + J R 2 + (z - L)' 



z + L + JR 2 + (z + LY 



(7) 



To obtain the simple form of equation ([7]) from the Legendre function Qq 
we used the auxiliar functions fi\ = z + L + Ri, fi2 = z — L + R2 and the 
identities [13] 



Ri 

z + L 



R 2 + £ 
2mi 



R 2 = 
z-L 



R 2 + ti 
_ A - R 2 

2^2 



(8) 
(9) 



A bar with maximum of density at the centre and vanishing density at both 
ends can be obtained by the superposition 



A02 — Ao — A2 
The corresponding potential reads 



^ m = ^E( R 2 +2L 2_ 2z 2 )ln 



3M 
4L 



(10) 



8L3 

3GM 
8L 3 



z-L+JR 2 + (z-L) 



z + L+\R 2 + {z + L) 



(L - 3z) JR 2 + (z + L) + (L + 3z) J R 2 + (z - L 



111) 



We will also consider another bar with density obtained by the superposition 

v2 



5 12 IbMz 2 
A024 — A + -A 2 - — A 4 — 



(12) 
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The density (|12p vanishes at the centre and at both ends of the bar, and has 
maxima at z/L = ±\/2/2. The associated potential can be expressed as 



$ 024 = (_ 3i? 4 + UR 2 Z 2 _ AR 2 L 2 + 8z 2 L 2 _ 



x In 



32L 5 

L+Jr 2 + {z-L) 2 



5GM , 2 _ 2L 2l 
32L 5 LV 



z + L + x /R 2 + (z + LY 

22zL 2 + 50z 3 - 6L 3 ) \/# 2 + (z - L) 



+ (55R 2 z - 9R 2 L + 26z 2 L + 22zL 2 - 50z 3 - 6L 3 ) J R 2 + (z + I? 2 



(13) 



The above potential-density pairs refer to infinitesimally thin bars, thus 
the potential is singular along the bar. For astrophysical applications (e.g. 
galactic bars) more realistic potentials should be free of singularities. A very 
simple way to 'soften' these potentials is by making a Plummer-like trans- 
formation R 2 — > R 2 + 6 2 , where b is a non- negative parameter [11] . With 
this procedure one obtains potential density-pairs that make a transition 
between infinitesimally thin bars (b = 0) and a Plummer sphere (6 S> L) [3]. 
Applying this transformation on the potentials 0, (fTTI) and (fT3|) . the cor- 
responding mass density distributions are calculated directly from Poisson 
equation in cylindrical coordinates, 

*RR + ^ + $,zz = 4vrGp. (14) 

The explicit expressions are given in Appendix [Aj The three mass densities 
are free from singularities and non-negative everywhere. For large values of 
R and z, the mass densities decay with (R 2 + z 2 ) _5//2 , as can be verified by 
an asymptotic expansion or simply by noting that in this limit the densities 
approach that of the Plummer sphere, which decays as (R 2 + z 2 ) -5 / 2 . Thus, 
in principle, they fall fast enough to put a clear cut-off and consider them 
as finite. In Figs [Ha) and (b) we show some isodensity contours of the 
dimensionless density po = po/(M/L 3 ), equation (j40j) . as functions of R/L 
and z/L for a 'softening parameter' b/L = 0.25 in Fig. Q^a) and b/L = 0.75 
in Fig. QJb). Figs a) and (b) and Figs[3]||a) and (b) display, respectively, 
isodensity contours of the other dimensionless barred densities (|42p and (|44p 
for the same values of the parameter b/L as in Figs QJa) and (b). The 
softened bars retain the same qualitative characteristics as the infinitesimally 
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(a) (b) 




012345 012345 
R/L R/L 



Figure 1: Isodensity contours of the dimensionless density po = po/iM/L 3 ), 
equation (|40|) . as functions of R/L and z/L for (a) b/L = 0.25 and (b) 
b/L = 0.75. 



thin ones, e.g. the isodensity curves in Figs[]Ja) and (b) are more elongated 
than those displayed in Figs[2]Ja) and (b), because the linear density of the 
thin bar (J6|) is less concentrated at its centre than the density of the thin 
bar (HD|). 



3 Equilibrium points and their stability 

An important aspect related to the morphology of barred galaxies is the 
study of motion of a test particle in the gravitational field of a uniform 
rotating bar. In this section we will discuss the equilibrium points and their 
stability for the motion in the field of the softened bars discussed in Section 
[5J For convenience, we place the bar along the rr-axis, and consider the 
motion on the xy plane. For this task, the potentials (f39l) . (j4"Tj) and (H3l) 
should be rewritten by replacing z — > x and R 2 — > y 1 . In the forthcoming 
discussion, we shall refer the potential-density pair (j39j) -(T10"1) as bar model 
1, the pair f)41f ) — (|42 |) as bar model 2 and the pair (|43p ~(|44 p as bar model 3. 

In a coordinate system attached to the bar that rotates with an (con- 
stant) angular velocity O, the equations of motion of a test particle are 

x-2ny = -^, (15) 

y + 2Qx = (16) 

dy 



5 



(a) 



(b) 
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Figure 2: Isodensity contours of the dimensionless density /?02 = 
P02/(M/L 3 ), equation as functions of and z/L for (a) 6/L = 0.25 
and (b) b/L = 0.75. 



(a) (b) 
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Figure 3: Isodensity contours of the dimensionless density /?024 = 
P02a/(M/L 3 ), equation (f¥i|) . as functions of and z/L for (a) 6/L = 0.25 
and (b) b/L = 0.75. 
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where dots indicate derivatives with respect to time, and & e ff. is the 'effec- 
tive' potential, 

®eff. = <S>bar ~ Y + ■ ( 17 ) 

At an equilibrium point, V<I> e //. = 0, and the resulting system of two alge- 
braic equations must be solved to obtain the equilibrium points (Lagrange 
points). Because of the symmetry of the models, the equilibrium points are 
symmetric with respect to the x- and y-axes. One Lagrange point is the ori- 
gin (0, 0), the pair on the x-axis will have coordinates (±xl, 0) and the pair 
on the y-axis will have coordinates (0, ±yi)- The stability of an equilibrium 
point is determined by the linearized equations of motion around it. The 
following conditions are necessary and sufficient for an equilibrium point be 
stable [3]: 

a/3 > 0, (18) 
- (a + /3 + 40 2 ) < 0, (19) 

(a + /3 + 40 2 ) 2 -4a/3 > 0, (20) 

where the second derivatives are evaluated at an equilibrium point. We 
shall calculate and analyse the stability of the Lagrange points for each of 
the three models of bars. 

3.1 Bar model 1 

At the origin the values of the second derivatives (|2ip are 

q = ° M , n -rf, f GU "» (22) 

One finds analytically that conditions (|19p and (j20|) are always satisfied. By 
condition (|18H . the origin will be unstable for angular velocities in the range 



GM GM 

< n < \hn^=^- ( 23 ) 



y(62 + L 2)3/2 ^vvwtt 1 

Fig. IHa) shows the stability diagram of the point (0, 0), as functions of b/L 
and of the dimensionless angular velocity f2 = ft/ W GM/L 3 . The unstable 
region grows as the bar becomes more elongated. 
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(a) 
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(b) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

b/L 

Figure 4: (a) Stability diagram of the equilibrium point (0, 0) for bar model 
1 as functions of b/L and Cl = Q/^GM/L 3 . (b) Stability diagram of the 
equilibrium point (0,ul/L) for bar model 1 as functions of b/L and Q, = 
n/y/GM/L 3 . 
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Figure 5: (a) Curves of equation (|26p for some values of b/L. (b) The 
stability of the equilibrium point(s) (xl/L, 0) for bar model 1 as function of 
b/L. 



On the y-axis the equilibrium point is given by the equation 



+ b2 ) Jvl + b 2 + L? 



(24) 



In this case the stability is better investigated by a graphical analysis of con- 
ditions (fl~8|) ~ ([20j) . Fig. LUb) displays the stability diagram of the Lagrange 
points on the y-axis as functions of b/L and Cl. For b/L ^ 0.75 the points are 
stable for all values of the angular velocity. For lower values of b/L there is 
an interval of Q where the equilibrium points are unstable and this interval 
becomes larger as b/L approaches zero. In this limit, the points are stable 
for an angular velocity less than 



U = 18 J 5 + 4V2\/ 2 + V2 ps 0.166. 

14 v v 

The equilibrium points on the x-axis are given by the equation 
GM 



(25) 



b 2 + {x L + L 



b 2 + (x L - LY 



- n 2 x L = (26) 
2L^P + (x L - Lf^jV + (x L + Lf 

For a given value of O the equilibrium point must be found by solving (I26p 
numerically. On the other hand, for a given value of xl one might calculate f2 
directly from (|26p . In Fig. EJa) we plot some curves of £1 as function of xl/L 
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for some values of the parameter b/L. It is seen that there may exist two 
equilibrium points for a given value of the angular velocity, which means 
two pairs of Lagrange points on the x-axis. This can happen for values 
of b/L ;$ 0.82. Fig. [5{b) shows the stability of the equilibrium point(s) 
(xl/L,0) as function of b/L. Comparing Figs a) and (b), we note that 
when two equilibrium points exist, the inner point is always stable, whereas 
the outer is unstable. When only one equilibrium point exists, it is always 
unstable. 



3.2 Bar model 2 



For this model of bar, the values of the second derivatives (f2Tj) at the origin 
read 



a 



3GM 
3GM 



In 



4L 3 



In 



Vb 2 + L 2 -L 
s Vb 2 + L 2 + L / 

'Vb 2 + L 2 -L S 
Vb 2 + L 2 + L 



3GM 



+ 



L 2 Vb 2 + L 2 
3GMVb 2 + L 2 



2b 2 L 2 



ti 2 



ti 2 



(27) 



The origin will be an unstable equilibrium point for angular velocities in the 
interval 



3GM 



In 




3GM 



L 2 Vb 2 + L 2 

^/b 2 + L 2 -L 
Vb 2 + L 2 + L 



1/2 



< n < 



3GM 
4.L 3 



+ 



3GM\/b 2 + L 2 

~Wl 2 



1/2 



(28) 



In Fig. Ufa) we display the stability diagram of the point (0, 0), as functions 
of b/L and of the angular velocity Cl. Also in this model the unstable region 
grows as the bar becomes more elongated. 

On the y-axis the equilibrium point is given by the equation 



3GM 
TP" 



In 



'y 2 L + b 2 + L 2 -L\ ZGMsJyl + b 2 + L 2 



ly 2 L + b 2 + L 2 + L / 



2L 2 {y 2 L + b 2 ) 



n 2 = 0. (29) 



The stability diagram of the Lagrange points on the y-axis is displayed in 
Fig. E|h). For b/L ^ 0.56 the points are stable for all values of the angular 
velocity. In the limit of infinitesimally thin bar, the points are stable for 
U < 0.25. 
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Figure 6: (a) Stability diagram of the equilibrium point (0, 0) for bar model 
2 as functions of b/L and Cl = Q/^GM/L 3 . (b) Stability diagram of the 
equilibrium point (0, Vl/L) for bar model 2 as functions of b/L and Q, = 
n/y/GM/L 3 . 
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The equilibrium point on the x-axis is calculated from 



3GM 
~2lJ 



In 



x L -L + \b 2 + (x L -LY 



x L + L+ \lb 2 + (x L + L 
3GM 



2x L L 3 



b 2 + (x L + Lf - Jb 2 + (x L - L) 2 



n 2 = o. (30) 



In this case there exists only one equilibrium point for a given value of the 
angular velocity, and we found that this point is always unstable. 



3.3 Bar model 3 

For bar model 3, the values of the second derivatives (|21|) at the origin read 



o 



IhGM (3b 2 + L 2 ) i / V6 2 + L 2 - L\ 15GM (3b 2 + 2L 2 ) ^ 2 



2L 5 



Vb 2 + L 2 + L 



L^b 2 + L 2 



(31) 



15GM (3b 2 + 2L 2 ) ( V6 2 + L 2 - L\ A5GMVb 2 + L 2 ^ 
p = In I . I — i I 



8L 5 



V6 2 + L 2 + L 



4L 4 



(32) 



The origin will be an unstable equilibrium point for angular velocities in the 
interval f^i < Q, < O2, where 



IbGM (3b 2 + L 2 ) ^ ( V6 2 + L' 2 — 15GM (3b 2 + 2L 2 ) 



2L 5 



Vb 2 + L 2 + L 



L^b 2 + L 2 



15GM (36 2 + 2L 2 ) / ^b 2 + L 2 - l\ 45GMV& 2 + L 2 
8L5 n \Vb 2 + L 2 + L J 4L4 



1/2 

1 

(33) 
1/2 

(34) 
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Furthermore, there is another region of instability for angular velocities 
greater then O = w A/B, where 



2\ 2 



.4 



2025G 2 M 2 (5b 2 + 2L 2 ) 



hi 2 



V6 2 + L 2 -L 



64L 10 

675G 2 M 2 (75b A + 85b 2 L 2 + 22L 4 



In 



\/6 2 + L 2 + L J 
f Vb 2 + L 2 -L 



16LV6 2 + £ 2 \ v v / 6 2 TX 2 + L / 

225G 2 M 2 (156 2 + 11L 2 ) 2 



15GM (96 2 + 2L 2 ) 

B = K — Mn 

L 5 



16L 8 {b 2 + L 2 ) 

' V6 2 + L 2 -L\ 30GM (9b 2 + 5L 2 ) 



Vb 2 + L 2 + L 



L^b 2 + L 2 



(35) 
(36) 



Figs 0(a) and (b) show the stability diagram of the point (0,0), as functions 
of b/L and the angular velocity £1. Fig. [T^b) gives an enlarged view of the 
second region of instability. 

On the y-axis the equilibrium point is given by the relation 



15GM (3y 2 L + 36 2 + 2L 2 
8£5 



In 



'yl + b 2 + L 2 -L 
!y 2 +b 2 + L 2 + L 



45GM^y 2 L + b 2 + L 2 
il 4 



Q 2 = 0. (37) 



The stability diagram of the Lagrange points on the y-axis is displayed in 
Fig. 0(c). For b/L ^ 0.90 the points are stable for all values of the angular 
velocity. In the limit of infinitesimally thin bar, the points are stable for 
Cl < 0.13. 

The equilibrium point on the x-axis is calculated from 



15GM (36 2 + L 2 - 2x 2 L ) 



2L 5 



In 



x L -L + Jb 2 + (x L -L) 



x L + L + Jb 2 + (x L + LY 



5GM (,~ r2 „ 2 , 
+ 2x^\ [ ' b+L llXL) 



b 2 + (x L + L) 2 -Jb 2 + (x L -L) 



-5Lxl 



b 2 + (x L + L) z + Jb 2 + (x L -L) 



n 2 



0. (38) 



13 



(a) (b) 




0.2 0.4 0.6 0.8 1 
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Figure 7: (a)-(b) Stability diagram of the equilibrium point (0, 0) for bar 
model 3 as functions of b/L and £l = fi/y GM/L 3 . (c) Stability diagram 
of the equilibrium point (0,yz,/L) for bar model 3 as functions of b/L and 

n = n/y/GM/L 3 . 
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Figure 8: (a) Curves of equation (|38p for some values of b/L. (b) The 
stability of the equilibrium point(s) (xl/L, 0) for bar model 3 as function 
of b/L. The dashed curve indicates the position of the static equilibrium 
point. 



Fig. [8][a) shows some curves of O as function of xl/L for some values of 
the parameter b/L. As happened with bar model 1, there may also exist 
two pairs of equilibrium points on the x-axis. This is possible for values 
of b/L ;$ 1.63. Fig. El^b) shows the stability of the equilibrium point(s) 
(xl/L,0) as function of b/L. Also here, when two equilibrium points exist, 
the inner point is always stable, whereas the outer is unstable. When only 
one equilibrium point exists, it is always unstable. 

From Fig. EJa) we note a particular feature of this model of bar: even 
without rotation (0 = 0), there is an equilibrium point along the x-axis for 
some values of the parameter b/L (for instance, x/L ~ 0.6 for b/L = 0.25). 
This static equilibrium point exists because the mass density of the bar is 
not concentrated at the origin (see Fig. [3]). In Fig. E^b) the dashed curve 
indicates the location of this point as function of b/L. A similar static 
equilibrium point was found in potential-density pairs for flat rings |14j . 

4 Discussion 

We presented analytical potential-density pairs for infinitesimally thin and 
softened bars constructed from an identity that relates multipolar solutions 
of the Einstein equations to Newtonian potentials of bars with densities pro- 
portional to Legendre polynomials. The main advantage of these models is 
that all potential density-pairs can be explicitly expressed in terms of ele- 
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mentary functions, and bars with a desired density profile can be constructed 
from the set of densities ([5]). 

As an application of the barred potentials, we calculated the equilibrium 
points for the motion of test particles in the gravitational field of three mod- 
els of rotating bars and analysed their stability. The results suggest some 
conclusions. The equilibrium point (0, 0) has the tendency to be more stable 
in bar model 2, and more unstable in bar model 3. The stability diagrams 
for the equilibrium points along the y-axis have the same qualitative be- 
haviour for the three models of bars. On the other hand, the properties of 
the equilibrium points along the rc-axis seem to be quite sensitive to the par- 
ticular model of bar used. In the case of bar model 2, the points are always 
unstable, whereas for bar models 1 and 3, there even exists the possibility 
of two pairs of equilibrium points, one being stable and the other unstable. 
It is known that the equilibrium points on the rr-axis of a homogeneous el- 
lipsoid are always unstable [U [6] . Prom our three models the bar model 2 
has the nearest shape of an ellipsoid, thus it is not surprising that it exhibits 
similar properties. It seems that barred mass distributions with less mass 
concentrated around the centre of the bar tend to stabilize the equilibrium 
points along the x-axis. Our results are in qualitative agreement with those 
obtained by Michalodimitrakis [6], who compared the stability properties 
of equilibrium points for a homogeneous ellipsoid and for a homogeneous 
parallelepiped. 
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A Potential- density pairs for softened bars 



The expressions for the potentials and mass densities p for the softened 
thin bar potentials ((7J), (fTT|) and (p~3|) are given by 



<I>r 



GM 



2L 



In 



z-L + K 2 \ 
z + L + Ux) 



Mb 2 



(z + L) [3 {R 2 + b 2 )+2{z + L) 2 



8irL(R 2 + b 2 ) 2 K\K\ 

-K\ (z-L) 3 (R 2 + b 2 )+2(z- L) 2 } , 

z-L + K 2 \ 
z + L + lli) 



02 



^L(R 2 + b 2 + 2L 2 -2z 2 )\n 



3GM 
' 8L 3 



(L - 3-z) ^1 + (L + 3z) ft 2 ] 



(39) 



(40) 



(41) 



3M6 2 



P02 



8ttL 3 {R 2 + b 2 yn 1 n 2 



(R 2 + 6 2 ) (^1 - K 2 ) z + ^1 (2 - L) 2 (s + L) 



-1Z 2 (z + L) 2 (z-L) 



024 



15GM 
32L 5 



(42) 



-3 (i? 2 + b 2 ) + 24 (i? 2 + 6 2 ) z 2 - 4 (# 2 + b 2 ) L 2 + 8z 2 L 2 

(S^rt) + w ( I- 55 < R2 + » 2 > 3 - 9 ^ + ft2 > i 

+26z 2 L - 22zL 2 + 50z 3 - 6L 3 ] 7£ 2 

+ [55 (R 2 + b 2 ) z-9 (R 2 + b 2 )L + 26z 2 L + 22zL 2 - 50z 3 - 6L 3 ] TZi) , 

(43) 



£>024 

{Tlx 



45Mb 2 
16vrL 5 
>2 , ,2\2 



In 



2 - L + ft 2 \ 



.z + L + ^J ' l6TTL 5 (R 2 + b 2 ) 2 TZ 1 TZ 2 
[R 2 + 6 2 ) z (5z + 3L) + (# 2 + b 2 ) (7z 3 - 5z 2 L - zL 2 + L 3 ) 



+ 



15M6 2 



+2z 2 (z + L) {z-Lf 



+ K 2 



{R 2 + b 2 ) 2 (-5z + 3L) + (R 2 + b 2 



-7z 3 - 5z 2 L + zL 2 + L 3 ) - 2z 2 (z -L){z + L) 2 } , (44) 
when- 7? , - \lR 2 + b 2 + (z + L) 2 and K 2 = J ' R 2 + b 2 + (z - L) 2 . 
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